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ABSTRACT 

A pathway to the formation of planetesimals, and eventually giant planets, may occur in concentra¬ 
tions of dust grains trapped in pressure maxima. Dramatic crescent-shaped dust concentrations have 
been seen in recent radio images at sub-mm wavelengths. These disk asymmetries could represent 
the initial phases of planet formation in the dust trap scenario, provided that grain sizes are spatially 
segregated. A testable prediction of azimuthal dust trapping is that progressively larger grains should 
be more sharply confined and furthermore the trapped grains should follow a distribution that is 
markedly different from the gas. However, gas tracers such as CO and the infrared emission from 
small grains are both very optically thick where the submm continuum originates, so observations 
have been unable to test the trapping predictions or to identify compact concentrations of larger 
grains required for planet formation by core-accretion. Here we report multifrequency observations of 
HD 142527, from 34 GHz to 700 GHz, that reveal a compact concentration of ^cm-sized grains, with 
a few Earth masses, embedded in a large-scale crescent of ^mm-sized particles. The emission peaks 
at wavelengths shorter than mm are optically thick and trace the temperature structure resulting 
from shadows cast by the inner regions. Given this temperature structure, we infer that the largest 
dust grains are concentrated in the 34 GHz clump. We conclude that dust trapping is efficient for 
approximately cm-sized grains and leads to enhanced concentrations, while the smaller grains largely 
reflect the gas distribution. 

Subject headings: Protoplanetary disks — Planet-disk interactions — Stars: individual: (HD 142527) 


1. INTRODUCTION 

Giant planet formation occurs in the first few million 
years following stellar birt h, while the parent pro toplan- 
etary disk is still gas-rich ( Zuckerman et al.|1995 ). How¬ 
ever, the classical debate on the formation mechanism, if 
envelope accretion o nto a rocky core, i.e. core-a ccretion 
( Pollack et al.|1996 ), or gravitational instability ([Kuiper 
/95l|), is stalled without observations. The pathway to 
giant planet formation determines the initial configura¬ 
tion and compositions of the product planetary systems. 
Current trends in theory contemplate a variety of forma¬ 
tion scenarios, notably second-generation core-accretion 
at large radii following the formation o f a giant closer-in 
(|Ayliffe et al.||2012| |Sdndor et al.||2011|). Models predict 
that young giant protoplanets carve a deep gap in the 
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dust componen t of protop l anet ary disk s, and a s hallow er 
gap in the g as (|Paardekooper & Mellema||2006 
et al. 2010|. The clearing of the protoplaneta 
2;hi 


Fouchet 


clearing of the protoplanetary gap is 
thought to underlie the class of ‘transition’ disks. A local 
pressure maximum develops at the outer ed ge of the gap, 


that can trap and pile-up the larger grains (Zhu & Stone 


2014), which would other wise rapidly migrate i nwards 
due to aerodynamic drag (|Weidenschilling 1977|) (as in 
the so-called ‘meter-size barrier 7 , which at 50-100 AU 
corresponds to mm-sized particles). Local pressure max¬ 
ima could also occur in l opsided disks with a stellar offset 
(|Mittal & Chiang][2015|. Whichever is origin, the devel¬ 
opment of a local pressure maximum may promote the 
pile-up and growth of dust grains, and the formation of 
large dust cl umps that could eventually collapse in plan¬ 
etary cores (|Ayliffe et al.||2012| |Lyra et al.||2009| |Sdndor 
et al. 2011). 


Recent observational progress brought by the Atacama 
Large Millimeter Array (ALMA) fits in the scenario of 
second-generation planet-formation at large stellocentric 
radii. The observational identification of so-called dust 
traps, in the form of concentrations of mm-size d dust 


(Casassus et al. 


2013 


van der Marel et al. 2013) in the 


outer regions of transition disks, suggests that azimuthal 
dust trapping may occur. An important prediction of 
this scenario is that pro gressively larger grains should be 


more sharply confined (Birnstiel et al. 


Lyra & Lin 


2013), possibly leading to the formation oFToulders and 
planetesimals. Another consequence is that the larger 
dust grains should have a distribution that is markedly 
different from the gas. Observed contrast ratios of 30 
to 100 between extrema in the sub-mm continuumm (at 
coarse angular resolutions) would seem too large to re- 
























































2 


Casassus et al. 


fleet an equally lopsided gas distribution. However, a 
theoretical question remains as to the impact of enhanced 
cooling in the dust trap on the physical conditions. In 
observational terms, the use of gas tracers such as CO 
or the infrared (IR) emission from small grains are both 
very optically thick where the ALMA thermal continuum 
originates, so that the previous observations ([Casassus 


et al.||2013 van der Marel et al.||2013 Perez et al.||2U14p 

could not test the trapping predictions! For instance, 

in IRS 48 th e outer ring is optica lly thick in the mid- 
IR emission (|Bruderer et al.||2014| their Sec. 4.1), as is 
CO(6-5) whose azimuthal structure is largely determined 
by foreground absorption. 

Since emission at a given wavelength is dominated by 
grains of a matching size, multi-wavelength radio obser¬ 
vations can confirm if the crescent-shaped continua seen 
in transition disks are indeed due to the dust trap phe¬ 
nomenon. In this work we compare new ALMA high- 
frequency radio continuum observations of HD 142527, 
obtained at 700 GHz (ALMA band 9; wave lengths of 
^0.43 mm, described in Appendix, Sec. A.l), with re¬ 
solved images at 34 GHz ( 8.8 mm) acquired at the Aus¬ 
tralia Telescope Compact Array (ATCA, see Sec. A. 2 ). 
In Section [ 2 ] we show that the optical depths are of or- 
der ^1 at 345 GHz, so that the sub-mm spectral trends 
observed with ALMA correspond to optical depths ef¬ 
fects rather than to dust trapping. In Sec. [3] we explain 
that the ATCA observations reveal a compact clump at 
34 GHz embedded in the sub-mm crescent, which is re¬ 
flected by an opacity spectral index corresponding to 
larger grains, as expected in the dust trapping scenario. 
In Sec. 0] we discuss the observational results in terms of 
radiative transfer predictions based on a parametric dust 
trap mode. Sec. [5] summarises our results. 

2 . AN OPTICALLY THICK SUB-MM CONTINUUM 

The ALMA band 9 continuum (Figure |T]) has a mor¬ 
phology that is intermediate bet ween the broken ring 
seen in the thermal IR at 18 /im (|Fujiwara et al. 2006; 
Verhoeff et al.| 2 011) and the 345 GH z crescent seen m 
AL1V1A band / (Casassus et al. 2013). The maximum 
radial width of the crescent, measured at half-maximum, 
is about 0.6", so is resolved with a beam size of 0.25". 

A traditional means to quantify the spectral variations 
as a function of position is the spectral index map a(x), 
with a power-law parameterisation of the specific inten¬ 
sity I v = l vo [y jvo)^^ . We find that the spectral in¬ 
dex between 345 GHz and 700 GHz, 0 : 345 , as well as the 
slope of the spectrum across the 0.13 GHz frequency lever 
within ALMA band 7, 0345 , both anti-correlate with the 
continuum crescent (see Sec. A. 6 ), meaning that the az¬ 
imuthal extension is greater at shorter wavelength. How¬ 
ever, while the band 7 and band 9 continua peak be¬ 
tween lOh and llh as a function of azimuth along the 
outer ring, the spectral index maps exhibit a common 
minimum at ^lh, i.e. at the northern ansa of the ring. 
This coincidence is indicative of a projection effect, as in 
limb-brightening: if column densities are greater at lh, 
the increased optical depths will result in flatter spectra. 

The spectral variations inferred from the multi fre¬ 
quency morphological trends can be cast into sky maps 
for the optical depth r(x) = r 0 (x) x (y/v 0 )PA x ) and 
line-of-sight temperature T s (x), as a function of angular 
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Figure 1 . ALMA band 9 observations of HD 142527, and 
comparison with synthetic predictions, x— and y— axes indi¬ 
cate angular offset in arcsec along right-ascension (RA) and 
declination (Dec) relative to the stellar position, at the origin 
of coordinates. The color scale shows a restored image of the 
ALMA band 9 at 700 GHz, with contours at 0.5 and 0.75 
times the peak intensity. The wedge indicates specific inten¬ 
sity in Jy beam -1 , with a beam of 0.21x0.17 arcsec (the beam 
ellipse is shown on the upper right corner). The inset shows 
emergent intensities predicted from the dust trap model (see 
Sec. El, including the impact on grain temperatures due to 
shadowing from a tilted inner disk. Contours are at 0.5 and 
0.75 times the peak, no smoothing has been applied. The side 
of the inset corresponds to 3.5 arcsec. 

position x. We fit the observations with the intensities 
eme rgent from a uniform slab (hereafter ‘grey-body’, see 
Sec. A.5 and Sec. 


C(^) = B v (T s (xj) [1 - exp(-r(f))]. 

The inferred temperature field approximates the opacity- 
weighted average temperature along the line of sight (see 
Sec. 4.1 for a discussion of biases). In synthesis, as ex¬ 
plained in Sec. [E| and illustrated in Fig. [8 the crescent 
reaches optical oepths close to 1 at 345 GHz, even in a 
rather coarse beam (which dilutes the signal). The opti¬ 
cal depth maximum is coincident with the location where 
the sub-mm spectral index (o™ and 0345 ) is minimum, 
thus in agreement with an interpretation of the spectral 
index variations in terms of optical depth effects, rather 
than genuine variations of the underlying dust grain pop¬ 
ulations. Another consequence of the high submm opti¬ 
cal depths is that the double-peaked crescent morphol¬ 
ogy, best seen in band 9, reflects structure in the tem¬ 
perature field rather than in the density field inside the 
dust trap. The local band 9 emission maxima, at 10.5h 
and 1.5h, are separated by a minimum in temperature at 
0.5h. 

3. ATCA/ALMA MULTI-FREQUENCY RESULTS 

3.1. A 34 GHz clump embedded in the sub-mm 
crescent 
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Resolved ATCA observations of HD 142527 in the opti¬ 
cally thin 34 GHz continuum reveal structure inside the 
sub-mm-opaque crescent, as shown in Fig. [2] (see also 
Fig. [3|. The lower levels at 34 GHz follow the Rayleigh- 
Jeans extension of the hotter dust seen in ALMA band 9. 
However, the peak signal stems from clumpy emission at 
^11.5h and a longer arc at lh to 2h. Both regions co¬ 
incide fairly closely with molecular decrements (Sec. |C|. 
Given the thermal noise in the field, the 11.5 h clump is 
significant at 10a, and it is robust against uv —filtering 
biases. To test for interferometer filtering artefacts, we 
ran Monte-Carlo simulations of ATCA observations on a 
deconvolved model image of the band 7 data (see Sec. 
thus bringing the multi-frequency data to an approxi¬ 
mately common Fourier basis. Fig. [2] shows that the 
ATCA clump stands out compared to the filtered ALMA 
maps. 

The peak 34 GHz intensity at 11.5h reaches 
147/iJy beam -1 , with a 0.64" x 0.30" beam, and has an 
equivalent brightness temperature of only 0.5 K, hence 
the signal is optically thin (or else it is beam-diluted). 
For the fiducial density of water, and for a single grain 
size of 1 cm with a circular cross section, at 20 K, we 
find that the observed flux density from the unresolved 
11.5h clump corresponds to >7.4 M 0 of dust. The in¬ 
equality accounts for the possibility of locking more mass 
in grains of different sizes that do not emit appreciably 
at 34 GHz. 

A compact signal coincident with the stellar position 
is conspicuous in the ATCA map. Its spectral index at 
ATCA frequencies is a^o = 1.0 ± 0.2, which extrapo¬ 
lates to the level of the faint circumstellar signal seen 
ALMA band 7 ( Casassus et al.||20l3 ; Fukagawa et al. 


m 


2013|). These spectral index values are consistent with 
tree-free emission associated with a stellar wind, or with 
stellar accretion. 


3.2. Evidence for dust trapping 

In another application of the grey-body diagnostics, 
we computed the optical depth and opacity index (/3s) 
maps corresponding to the ATCA and ALMA datasets 
shown in Fi g. |2} i. The resulting optical depth map 
shown in Fig. |2p confirms the levels at 345 GHz inferred 
from the ALMA band 7 and band 9 comparison. The 
ATCA clump at 11.5h translates to a minimum in f3s(x) 
(Fig#). This is an indication that dust grains in this 
clump are larger than in the rest of the crescent. As 
summarised in Fig. [9j the variations in /3s are significant 
at ^7 a. 

4. PREDICTIONS BASED ON A SYNTHETIC 
DUST TRAP MODEL 

4.1. Biases of the grey-body diagnostics 

The grey-body line-of-sight diagnostic is based on 2- 
dimensional sky maps for /3s(x) and Ts(x ), while the 
actual physical conditions are 3-dimensional. Since the 
crescent is optically thick in band 9, the temperature is 
mainly set by the higher frequencies, and probably over¬ 
estimates the temperature in regions closer to the mid¬ 
plane. Hence the actual optical depths could be higher 
than inferred from grey-body diagnostics. In order to 
assess biases and test the dust trap interpretation, we 
applied the same observational diagnostic, based on the 


grey-body fits, to synthetic radiative transfer predictions 
of a parametrised dust trap model meant to approximate 
the phenomenon observed in HD 142527. The synthetic 
model was implemented in RADMC3D (Dullemond et al. 
2015|) (Methods), and centers the dust trap at 11.5h. 
At the native resolution of the radiative transfer, the 
grey-body optical depth map reproduces the input op¬ 
tical depths very closely. However, both optical depth 
and temperatures are biased downwards by ^50% af¬ 
ter filtering by the ATCA response. The location of the 
minimum in the grey-body opacity index, / 3s(x ), coin¬ 
cides with the input dust trap, which is fairly optically 
thick in band 7. However the optical depths in the syn¬ 
thetic trap are not high enough to result in the observed 
deep molecular decrements - their properties could be 
reproduced by prescribing gas temperatures cooler than 
the dust (Sec. O. In particular the existence of decre¬ 
ments in CO(2Mj (Perez et al. 2015), so at frequencies 
with lower optical depths, suggests that there is ample 
room for other effects. The decrements may also corre¬ 
spond to structural features of the disk, i.e. shadowing 
in lower scale heights under the 34 GHz clumps, due to 
midplanes cooler than in the rest of the outer disk at the 
same stellocentric radius. A mechanism for such a local 
‘pinching’ of the outer disk may be found in the dust 
trap scenario. In radiative equilibrium the larger grains 
tend to be cooler, so that hydrostatic scale heights will be 
lower under the dust trap. The shadowed CO is cooler, 
or possibly frozen out and condensed onto grain surfaces. 

4.2. Temperature decrement due to shadowing on the 
outer-ring 

The temperature structure inferred from the band 9 
observations are readily explained by the model, which 
includes a tilted inner disk within 20 AU, as implied by 
the sh adows in scat tered light cast on the outer disk 
(Marino et al.||2015|). Fig. fl] shows that the predicted 
emergent intensities in band""9 are double-peaked, much 
like the observations. The warped structure results in 
shadows projected on the outer disk where the stellar 
heating is blocked by the inner disk at PAs of 0.5h and 
6.5h. The similarity with the observations (Fig. [8| is re¬ 
markable considering the idealizations of the synthetic 
disk, which assumes a circular cavity, power-laws for the 
grain size distributions, and steady-state dust trapping 
prescriptions. Another idealisation of the model is steady 
state passive heating. 

The cooling timescale is tc = U/(gbT a ), where U is 
the internal energy per unit area, and gb is the Stefan- 
Boltzmann cost ant. The condition that r c must be 
shorter than the crossing t ime under shadows that cover 
^20°at 140 AU ( Marino et al.| 2015, and Fig. 11), is sat¬ 
isfied for surface densities lower than 50 g cnT 2 , which 
is indeed the maximum surface density in the model. 
It is thus reasonable to expect that the double-peaked 
structure of the continuum emission is determined by the 
temperature field, and that the mm-grains are otherwise 
uniformly spread along the crescent. The modulation 
in temperature along the crescent, caused by the inner 
warp, impacts on the thermodynamics of the disk. This 
modulation imparts a periodic forcing to the outer disk, 
whose dynamical consequences should be investigated, 
and that could perhaps explain the observed spiral pat¬ 
tern. 
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Figure 2. Multi-frequency continuum highlighting larger grains at 11. 5h. The data have been filtered for the ATCA response, a: 
RGB image with the ATCA image in red (with corresponding contours at 0.4 and 0.7 times the peak intensity at 147/xJy beam -1 ), 
ALMA band 7 in green, and ALMA band 9 in blue (see Fig. [6] for a panel with the same shown images separately). The inset 
shows emergent intensities predicted from the dust trap monel, after filtering for the ATCA response, with 34 GHz in red 
contours at 0.5 and 0.75 times the peak. The side of the inset corresponds to 3.5 arcsec. b: Observed grey-body temperature 
map. c: optical depths at 345 GHz. d: emissivity index map /3 s (x), with ATCA contours in red. The /3 S minimum at 11. 5h 
indicates the presence of larger grains. 


4.3. Very lopsided disks? 

The spectral trends in HD 142527 at ALMA frequen¬ 
cies can thus be accounted for with optical depth effects, 
such that no trapping is required for grains up to 1 mm in 
size. In other words, the large sub-mm crescent (Casas 
sus et al.||2013|) mostly reflects the gas background, with 
relatively inefficient trapping, so that the observed con¬ 
trast ratio of ^30 is accounted for with a contrast of 20 
in the gas. Marked gaseous asymmetries with contrasts 
of order 10 are seen in MHD sim ulati ons of dynamical 
clearing by planet formation (Zhu & Stone 2014), and 
could perhaps reach higher values in models incorporat¬ 
ing dea d-zones with discontinuous turbulence prescrip¬ 
tions (Regaly et al. 2012), or with enhanced cooling in¬ 
side the"dust trap via dust-gas coupling. Thus we pro¬ 
pose that other sub-mm asymmetries could partly reflect 
the gas, perhaps even with contrasts of order ^100 such 
as in IRS 48, and most likely in the more mode rate c on¬ 
trasts , such as in HD 135344B and SR 21 (Perez et al. 
2014). 


5. CONCLUSION 


In summary, the multi-frequency observations of 
HD 142527 are consistent with the dust trap scenario. 
New ALMA observations have shown that the mm- 
grains are so abundant in the crescent that efforts to ob¬ 
serve grain growth at sub-mm wavelength are thwarted 
by optical thickness. Optical depth effects account for 
the spectral trends at ALMA frequencies, while shad¬ 
owing of UV-optical radiation by the central warp ex¬ 
plains the temperature structure of the crescent. How¬ 
ever, at transparent cm-wavelengths a compact clump of 
approximately cm-sized grains is embedded within the 
broader crescent of smaller mm-grains (the exact grain 


sizes depends on the size distributions). These phenom¬ 
ena are expected in the context of the aerodynamic cou¬ 
pling of dust grains with pressure gradients, given the 
physical conditions in the passively heated outer disk of 
HD 142527. Prospects are good for witnessing planetes- 
imal formation by dust trapping with radio observations 
of HD 142527 at cm-wavelength. 
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APPENDIX 
A. OBSERVATIONS 

HD 142527, at ^140 pc, is a promising laboratory to 
observe on-going gas giant planet formation, and con- 
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firm the existence of the dust trap phenomenon. This 
transition disk has a close to face-on orientation. In- 
fr ared^observatio ns found an inner disk, ^10 AU (Fuji- 


larly large gap (Fiika£ 

s;awa et al. 2006) 

disk (Casassus et al. 

2012 

| beyond 1 


tion. The cavity is dev oid of bodies larger than Mj up 


(Casassus et al. 2013), but an accreting companion has 


been reported at ^10 AU ( Biller et al.||2012 Close et al. 


2014 Rodigas et al.||2014) (with a binary mass ratio of 
0.01to 0.1). The bulk of th e mass lies in the outer disk, 


whos e crescent morphology ( Qhashi|2008 Casassus et al. 


emise 

Peres 

arms 

don dCasassus et al. 2013 van der Plas et al. 2014 

5 et al. 2015|. A complex array c 
sprouts away from the outer disk ( 

>f trailing spiral 

Fukagawa et al. 

2006; Casassus et al. 2012; 

Rameau et al 

. 2012; Canovas 

et al. 

20131 Avenhaus et al. 

2014; Christiaens et al. 2014). 

In 

this Section we provid 

Le details on the datasets that 


sustain our multi-frequency analy sis. Sec. ED presents 
the new ALMA obser vation s, Sec. |A.2| presents the new 
ATCA data and Sec. |A.3| su mmarizes previously pub¬ 
lished ALMA data. Sec. | A. 4| b riefly describ es ou r image 
synthesis procedure, while Sec. |A.5| and Sec. [AT] describe 
the alignment of the multi-frequency A LMA data and as¬ 
sociated spectral index trends and Sec. |A.7| describes the 


alignment of the ATCA and ALMA datasets. 

A.l. ALMA band 9 observations 


ALMA Band 9 observations of HD 142527 for program 
JAO.ALMA#2011.0.00465.S were carried out in 2012, in 
the nights of June 3, from 02:52:53 to 04:39:10 UT, and 
on June 4, from 00:43:22 to 02:42:26 UT. The precip- 
itable water vapor (PWV) in the atmosphere was stable 
between 0.54 mm and 0.62 mm on June 3, with a me¬ 
dian value at zenith of 0.581 mm. Conditions on June 4 
were similar, with a median PWV at zenith of 0.597 mm. 
The ALMA correlator provided 1875 MHz bandwidth in 
four different spectral windows at 488.28125 kHz resolu¬ 
tion (or 211.67 m s -1 at 691.47308 GHz) per channel. 
Each spectral window was positioned in order to target 
the CO(6-5) transition at 691.47308 GHz, HCO+(8-7), 
H 13 CO + (8-7), and HCN(8-7). Only CO(6-5) was de¬ 
tected. This line is the topic of a separate article. On 
June 3, the measured system temperatures ranged from 
~ 1400 zb 500 K in the spectral window covering CO(6- 
5) to ^ 2000 ± 700 K in the spectral window cover¬ 
ing HCO + . System temperature values 100 to 200 K 
higher were recorded on June 4. The number of 12 m 
antennas available at the time of the observation was 20, 
although one antenna reported very large system tem¬ 
peratures (DV12), and was flagged during data reduc¬ 
tion. Excluding calibration overheads, the total time 
on source for HD 142527 was 43 min on June 3, and 
49 min on June 4. The primary flux calibrator was Ti¬ 
tan, while 3C279 was used as bandpass calibrator. We 
used J1517-243 as phase calibrator, for which we measure 
a flux density of 0.70±0.14 Jy at 702 GHz, where the 
uncertainty accounts for systematic calibration errors. 
This measurement is consistent with historical values for 
this source. The radiometer phase correction errors were 
large enough to decorrelate the signal. We used the self¬ 
calibration algorithm to determine improved antenna- 


based phases that were consistent with the continuum 
image of HD 142527. 


A.2. ATCA data 

We used ATCA to observe HD 142527 in several ar¬ 
ray configurations in May 2010 (6A), June 2010 (6C), 
July 2013 (6A) and August 2013 (H168), with maximum 
baselines extending to 6 km. T he Compact A rray Broad¬ 
band Backend, or CABB (Wilson et al.| 2011|), provides 
two sidebands with bandwidths of 2 GHz in 2024 chan¬ 
nels each. We centred each sidebands at 19/21, 23/25, 
33/35 and 41/43 GHz (in this work we focus on the 33/35 
dataset). Complex gains were derived from observations 
of the quasar 1600-44, separated by about 2.8 degrees on 
the sky from HD 142527. The gain calibrator was typi¬ 
cally observed every 5 to 15 min for between 1 and 3 min 
duration, dependent on atmospheric conditions. Point¬ 
ing checks were also made on the quasar every 60-90 min. 
The bandpass response was determined from 15 minute 
observations of the quasar 1253-055 (3C 279). The ab¬ 
solute flux calibration was performed using ATCA’s pri¬ 
mary cm-band flux calibrator, the quasar 1934-638. All 
subsequent data reduction aud calib ration was performed 
with the Miriad software (|Sault et al. 19951). The CABB 
setup at 33/35GHz was used in 4 runs with different ar¬ 
ray configurations, including HI68, which allows us to 
estimate that the amount of flux loss in the extended 
configurations is negligible. The flux densities obtained 
in H168, 1.28±0.04/1.48±0.04, in Jy, for each sideband 
centered at 33/35GHz, is indistinguishable from the val¬ 
ues of 1.23±0.05/1.49±0.06, 1.34±0.11/1.48±0.11 and 
1.25±0.13/1.24±0.09 obtained in 6C, 6C and 6A, respec¬ 
tively. 

The ATCA data with the best dynamic range was ob¬ 
tained in Ka, i.e. at 33 GHz / 35 GHz. For efficiency 
we averaged each CABB 2 GHz spectral window into 
8 channels. We then brought all of the 6 km visibility 
data in Ka to a common epoch, which we chose as July 
2012, by correcting for the source proper motion with a 
translation phase (using the Miriad task uvedit). The 
6 km data in Ka w ere then concatenated with the task 
concat from CASA (McMulli n et al.| [2007) version 4.2.2. 
A standard clean restoration revealed a taint modulation 
at low spatial frequencies in the field, atypical of image 
synthesis artifacts and reminiscent of phase calibration 
errors. We therefore applied selfcalibration of the visi¬ 
bility phases, using the clean model and the CASA tasks 
gaincal and applycal. We used a solution interval of 
30 min and combined both spectral windows for sensi¬ 
tivity. The clean restoration of the corrected visibili¬ 
ties results in essentially the same source morphology as 
the uncorrected dataset, but with less systematics in the 
field. The impact of self-calibration can be assessed by 
inspection of Fig. [3j where we show the restored i mages 
obtained with our uvmem (Casassus et al. 2006[ 2013) 
tool before and after self-cali bration ( more information 
on this image synthesis is given in Sec. A. 4| below). 


A.3. Previous ALMA band 7 observations 

As explained in Fig. [4] and Sec. |A.5[ the calibrated 
dataset ALMA#2011.0.00465.S is affected by faint out¬ 
liers. We therefore used the ALMA Band 7 ob¬ 
tained by Fukagawa et al. 


(Fukagawa et al. 2013) 








































































Figure 3. Comparison of the restored ATCA images and the impact of self-calibration. The wedges indicate specific intensity in 
mjy beam -1 , a) no self-calibration ( uvmem restoration with a beam of 0.570 x 0.29 6 arcsec 2 ). b) self-ca librated image with uvmem 
restoration (with a beam of 0.639 X 0.302 arcsec 2 ). c) self-calibrated image with CASA (McMullin et al~] 2007) clean restoration (with a 
beam of 0.644 x 0.304 arcsec 2 , CASA version 4.2.2). 


ALMA#2011.0.00318. S, keeping only the last scheduling 
block, observed on 10-Aug-2012. The reason for not com¬ 
bining the whole dataset is because the dynamic range 
under the crescent is limited by image synthesis, rather 
than thermal noise, and because pointing uncertainties 
would degrade the beam if co mbine d without a tedious 
alignment procedure (see Sec. A.5). For instance, we 
measured an offset of 0.065 arcsec between this schedul¬ 
ing block and the single scheduling block in band 7 from 
ALMA#2011.0.00465. S, comparing exclusively the com¬ 
mon spectral window at 342.8 GHz. These offsets are 
within the nominal pointing accuracy of ~0.1 arcsec. Al¬ 
though small, such offsets have appreciable consequences 
when deriving spectral index maps. 

In summary, we used the dataset acquir ed by Fuk- 
agawa et al. (Fukagawa et al. 2013), program 
J AO. ALMA#2011.0.00318. S but tied the astrometric 
point ing to the data prev iously presented by Casassus et 
al. flCasassus et al.l|2013| ), JAO.ALMA#2011.0.00465. S. 
Furthe r det ails onth e astrometric registrations are given 
in Sec. IA.5I and IA.7I 


A.4. Image synthesis 

For image synthesis we use a non-parametric least- 
squares modeling technique (|Casassus et al. 2006) with 
a regularizing entropy term (i.e. as In the"family of max¬ 
imum entropy methods, MEM). We call this tool uvmem. 
These deconvolved model images are ‘restored’ by con¬ 
volution with the clean beam and by adding the dirty 
map of the visibility residuals calculated using the CASA 
package, for a chosen weighting scheme. In this work we 
used natural weights. 

The clean algorithm is an efficient and widely 
accepted technique to produce images from radio¬ 
interferometer data. To produce consensus images that 
can be compared with our own uvmem tools we used stan¬ 
dard CASA clean, i.e. Cotton-Schwabb. A comparison 
between uvmem and clean is shown in Fig. [3] 

We assume that the relative values of the visibility 


weights in the ALMA data are accurate. However, the 
absolute values do not match the scatter in the visibility 
samples recorded at each integration. In order to ex¬ 
tract statistics and perform least squares fits, we there¬ 
fore scaled the weights so that they correspond to 1/cr 2 , 
where a is the root mean square dispersion of all the 
visibility samples recorded for a given baseline and chan¬ 
nel. This procedure was carried out using a tool we call 
uvreweight. A single scale factor was applied to the 
weights of the whole dataset, given by the median val¬ 
ues of the ratios between observed scatter and tabulated 
weights. Thus we respected the relative values of the 
weights as provided in the calibrated dataset. We used 
median averages for a robust estimate of the observed 
scatters, after correction by the median absolute devi¬ 
ation. In this uvreweight-processed data the absolute 
value of the weights can be used to extract statistics. 


A.5. Alignment of the multi-frequency ALMA data 

We aligned the multi-frequency interferometer data by 
cross-correlating the visibilities, as illustrated in Fig. [4j 
We first approximately bring two independent visibility 
datasets to a common ^-coverage in the following way. 
We select as reference the dataset with the more compact 
uv -coverage (usually the lower frequency), and produce a 
deconvolved model using uvmem of the visibility dataset 
with the more extended uv— coverage. We then com¬ 
pute synthetic interferometer data corresponding to the 
more compact uv —coverage using our tool uvsim. For 
consistency the procedure is also applied to the refer¬ 
ence. Having produced comparable visibility datasets, 
we regularly sample shifts in the world coordinate sys¬ 
tems (WCS) of the two frequencies in a uniform grid of 
positional offsets (assuming that the parallactic cal¬ 
ibration is perfect in both cases).For each WCS offset Si 
we computed shifted band 9 visibilities 


r 6 s 9 hlft (ufc) = V b9 (u k )ex p 


z27r(<5j • Uk) 
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Figure 4. Band 7 vs. Band 9 visibility scatter plot, from 
dataset 2011.0.00465. S. Top: cross-correlation before out¬ 
lier filtering. Middle: after outlier filtering, but before as¬ 
trometric WCS correction. Lowe r plot: after the alignment 
procedure described in Sec. |A.5| In this case, the applied 
WCS offset is (-0.063", -O.Ol^Talong R.A. and Dec. We see 
here that this band 7 dataset ALMA^2011.0.00465.S is af¬ 
fected by faint outliers, t his is why we used the data from 
ALMA#2011.0.00318.S (Fukagawa et al.||2013l) . Since the 
dynamic range under the crescent is not limited by thermal 
noise, we kept only the last scheduling block, observed on 10- 
Aug-2012, to minimize pointing offsets. For consistency we 
tied the astrometric pointing to 2011.0.00465. S. 


We then calculated a linear regression between real and 
imaginary parts of visibilities, used as independent data, 
and computed the L 2 distance (squared norm) between 

model and data. We form an image X 2 ({^i}) of the L 2 
distances for each 5{. The best WCS shift is given by the 
centroid 5 0 of an elliptical fit to X 2 ({^})- Although we 
measure the shifts S Q in the uv— plane, the actual appli¬ 
cation of S 0 to bring images to a common WCS system 
is carried out after image restoration, so in the image 
plane. 


A.6. Spectral index maps 

Thanks to the positional alignment procedure, we can 
build a Band 7/Band 9 spectral index map by bringing 
both datasets to an app roximately common nu-coverage, 
as outlined in Sec. |A.5| We then compute restored im¬ 
ages of the uvmem models. The resulting Band 7/Band 9 
spectral index map is shown in Fig. [5) We caution that 
the absolute values of the spectral index map a is depen¬ 
dent on the chosen absolute flux calibration and SED, 
as explained in Sec. [D| However, the relative variations 
in a are significant, and independent of the chosen scale. 
Fig. [5] corresponds to the flux scale labelled CISscale (see 
Sec.]D|). 

In order to build a spectral index map based on the 
intra-band frequency lever of the ALMA band 7 data, 
we selected the lowest frequency spectral window as ref¬ 
erence 'uu-coverage. For instance in the spectral fits to 
the intra-band spectral windows (SPWs) of the ALMA 
band 7 data ALMA#2011.0.00318.S, this corresponds to 
329.3 GHz. We then bring all 4 spectral windows in 
band 7 to a comparable uv- c overage, by following the 
procedure outlined in Sec. |A.5| (so the reference SPW was 
also processed in the same way as the other 3 SPWs). 

Once all four spectral windows are brought to a com¬ 
parable uv— coverage, we extract a spectral index map by 
fitting the restored specific intensity maps with a power- 
law. There are therefore two free parameters per pixel, 
the intensity at the reference frequency, and the power- 
law exponent. Fig. [5] shows the resulting fit, along with 
the 1 a error on the spectral index. We caution that the 
absolute values of this intra-band spectral index map de¬ 
pends on the choice of SED scale, as explained in Sec.[D 
However, the morphological trends are independent of 
the total flux density scale. 


A.7. Alignment of the ATCA and ALMA data 

The offset between the ATCA and the ALMA Band 7 
data from ALMA#2011.0.00465.S was within our mea¬ 
surement errors. Since in the ATCA data we have 
two IFs centered on 33 GHz and 35 GHz, we can test 
the alignm ent by cross-correlating each IF separately, 
as in Sec. [A3J and also by requiring consistency be¬ 
tween real and imaginary parts. It turned out that 
the ATCA/Band 7 offset is systematically smaller than 
0.15", and changes sign. We concluded that the ATCA 
and ALMA pointing are consistent, but that the mor¬ 
phologies are too different for a more refined alignment 
with our cross-correlation scheme. We also attempted 
the standard strategy of identifying a common anchor 
point, in this case the stellar position. We document 
this exercise as it bears on the morphology of the sig¬ 
nal from the star’s vicinity in ALMA band 7. The cen¬ 
tral free-free signal, at the limit of dynamic range in 
band 7, is difficult to use as astrometric reference as 
its morphology is sensitive to self-calibration, and va ries 
from a compact inner disk reported by Fukag awa et al. 
(2013 ) to a gap-crossing filament in our previous analy¬ 
sis flCasassus et al7||2Q13 ). However, the continuum fila¬ 
ment’"weTmage" in our self-calibrated band 7 data, after 
filtering for outliers using our tool uvreweight, culmi¬ 
nates in a local maxim um near the cent r al sta r, now in 
closer agreement with Fukagawa et al. (2013). Given 
the uncertainties affecting the very taint stellar signal in 
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Figure 5. Spectral index map from the aligned images of the ALMA Band 7 and Band 9 datasets, a) spectral index from the 
comparison between band 7 and band 9, a™. The single contour shows the Band 7 continuum traced at half-ma ximum, b-d) 
Power-law fits to the intra-band 7 ALMA data, using dataset ALMA#2011.0.00318.S (Fukagawa et al.||2013|. b) spectral 
index map a(x). c) best-fit amplitude at the reference frequency of 345 GHz, in Jy/sr. d) thermal uncertainty on the spectral 
index map. 


band 7, the default astrometries for ATCA and ALMA 
band 7 ALMA#2011.0.00465.S appear to be consistent. 
We note, however, that the 11.5h emission appears to 
extend outwards in radius, by about 0.3 db 0.1" com¬ 
pared to ALMA band 7, while the lh clump does not. 
But given the above pointing uncertainties, this effect is 
not significant enough to warrant further consideration 
and modeling. We brought the multi-frequency data of 
HD 142527 to a common uv —coverage following the same 
procedure as for the multi-frequency ALMA data. The 
main difference is that here we used the ATCA data as 
reference. The resulting restored images are shown in 
Fig. [2ji, where the origin of coordinates corresponds to 
the location of HD 142527 at the Jun 2012 epoch. The 
central point source at 34 GHz is satisfactorily aligned 
with the origin. 

B. MONTE-CARLO TESTS FOR THE 
MULTI-FREQUENCY MORPHOLOGICAL 
DIFFERENCES 

Images obtained with pre-ALMA radio-interferometry 
are known to be affected by w-coverage artifacts, espe¬ 
cially with ATCA. Despite the efforts invested in cov¬ 
ering the i^v-plane in several ATCA runs, with different 
array configurations, it is important to test for such uv- 
coverage artifacts, that could for instance inject the 11.5h 
local maximum. We first produce deconvolved model 
images of the ALMA band 7 and band 9 visibility data 
using our tool uvmem, and then simulate ATCA observa¬ 
tions on these model images, following exactly the same 
'Mf-coverage as in the ATCA 34GHz concatenated data 
using our tool uvsim. The simulated visibilities are then 
scaled by a linear regression (forced to cross the origin), 
so that the simulated and ATCA visibilities are com¬ 
parable in magnitude. Following this scaling, we inject 
thermal noise given by the observed visibilities weights 
(having first applied our tool uvreweight to replace the 
ATCA visibility weights by the observed scatter). We 
repeat this procedure 200 times before taking statistics. 

Fig. [6] summarizes the result of these Monte-Carlo sim¬ 
ulations. We see that the signal is more extended in the 


ALMA bands than it is in ATCA, so that the visibility 
scaling we applied results in lower signal-to-noise ratio 
than in the actual ATCA observations. In this sense our 
MC simulations are conservative, we are injecting more 
noise relative to the signal (by about 40% in band 7). 
The final rms dispersions are up to ~1 /7 of the differ¬ 
ences between the ATCA and ALMA images - in other 
words, the differences between ATCA the ALMA images 
range from — 3cr to +6.2a. The reduced % 2 taken over 
pixels where both ATCA and ALMA are above 1/3 of 
the peak, is 4.2 for band 7, and 4.9 for band 9, for 12 
independent data points (the number of ATCA beams 
enclosed in the 0.3 times maximum contour) - the im¬ 
ages are thus different with a very high confidence level. 
In particular, under the 11.5h clump (so at the ATCA 
maximum), the difference with the band 7 average is 5.9 
times the band 7 rms. 

An additional test was performed by subtracting the 
synthetic band 7 ATCA visibilities from the 34 GHz 
ATCA data, after appropriate scaling by linear regres¬ 
sion. The stellar emission stands out in the residuals, as 
well as the llh ATCA clump. These significant resid¬ 
uals, both positive and negative, point to genuine mor¬ 
phological differences between 34 GHz and 345 GHz. In 
particular, the peak residual corresponding to the llh 
clump is at 50 /iJy beam -1 , while the thermal noise is 
at 10 fi Jy beam -1 . 


C. COINCIDENCE BETWEEN ATCA MAXIMA 
AND MOLECULAR LINE DECREMENTS 

The 34 GHz clumps coincide with the molecular-line 
decrements in HCO + (4 — 3) a nd in CO(3 — 2) that have 
been previously documented (Casassus et al. 2013). It 
has been suggested that an optically thick continuum 
may somehow be responsible for the correlation of th e 
12 CO decrements with the 345 GH z crescent (Casassus 
|et al 12013 van der Plas et al.|2014 ). Indeed, as shown in 
Fig.[?J the synthetic dust trap model documented below 
readily reproduce such decrements with a dust to gas 
mass ratio of 94, a standard CO abundance of 10 -4 , and 
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Figure 6. Statistics of the multi-frequency morphological differences inferred from Monte-Carlo (MC) simulations. All images 
have been restored from uvmem models, a): restored ATCA image at 34 GHz. b): average of ATCA observations with 200 
realizations of noise on the band 7 data, c): average of the MC simulations for band 9 d)-g): different realizations of noise, 
h): rms scatter of the band 9 simulations, i): rms scatter of the band 7 simulations. 


the following temperature profile: 

T = 70 K, if r < 30 AU, or (Cl) 

T = 70 K(r/(30 AU) -0 ' 5 if r > 30 AU. 

However, setting the gas temperature to a representative 
dust species in the outer disk (we chose the average popu¬ 
lation of amorphous carbon grains) attenuates the depth 
of the decrements - they are restricted to the outline of 


the underlying optical depth field. The decrements are 
reproduced by the above prescription for the gas temper¬ 
ature because the gas temperatures in the crescent, at 
r = 140 au, are ^30 K, and cooler than the continuum 
temperatures in the model, which reach 50 to 60 K. Thus 
in this case the gas acts like a cold foreground on a hot¬ 
ter continuum background, which after continuum sub¬ 
traction appears like foreground absorption. In fact, CO 
at higher altitudes could be hotter than the background 
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Figure 7. Frequency-integrated image of the 12 CO (3-2) emis¬ 
sion calculated with RADMC3D on our preferred model for 
the dust trap of HD 142527. x— and y- axis indicate offset 
from the stellar position along R.A. and Dec., in arcsec. 


continuum, so could still stand out after the subtraction 
of a colder optically thick continuum. Yet the gas tem¬ 
perature prescription above is inspired by the molecular 
line data, and higher temperatures can be ruled out as 
they greatly overestimate the integrated CO line fluxes. 


D. ABSOLUTE FLUX SCALE 


The inferred grey-body diagnostics and spectral index 
maps depend on the absolute calibration of the input 
multi-wavelength data. There are however substantial 
uncertainties on the observed SED. For instance, the 
flux densities reported for di fferent ALMA programs in 


band 7 vary by up to 50%, (Casassus et al. 


agawa et al. 2013). The band 9 flux densities are af- 


2013; Fuk- 


tected by similarly l arge uncertainties. The H erschel 
SPIRE flux density (van der Wiel et al.| 2014) at the 
same frequency is 13.7 Jy. We investigated two plausible 
scenarios for the observed SED. The pectral trends re¬ 
ported here are derived from the first of these flux scales, 
which we label CISscale , where we adopted a reference 
band 7 flux density of 3.5 Jy at 345 GHz, and took the 
SPIRE flux density for band 9. The comparison between 
the SPIRE flux density and the band 7 levels implies 
a spectral index ~ 1.9. O n the o th er han d, the 
band 6 ALMA data at 250 GHz ( Perez et al.| 2015) yields 


a 


2.7. We thus assumed that a plausible value for 
the intra-band 7 spectral index is a^ = 2.2. 

The synthetic dust trap, documented below, fits the 
SED within the error bars. So, for the purpose of com¬ 
paring the observations with the synthetic dust trap pre¬ 
dictions, we also considered a second flux sale, which we 
labelled SEDscale. This scale has a rather high band 9 
flux density of 15.7 Jy at 700 GHz, and is rather low 
in band 7, with 2.17 Jy at 342.9 GHz, while it comes 
very close to the reported flux density at 34 GHz, with 
1.34 mJy. All morphological trends are similar in both 
flux scales, but the highest optical depths in the SEDscale 
are lower - dropping by about 2.5 in CISscale to 0.7 in 
SEDscale (at 345 GHz). The temperature map appears 
to be more irregular in SEDscale , which suggests that 
CISscale may be more realistic. 


E. LINE OF SIGHT TEMPERATURES AND 
OPTICAL DEPTHS 

The grey-body uniform-slab models are used as diag¬ 
nostics for the optical depth, and are meant to charac¬ 
terize the emergent spectrum at a given line of sight with 
three parameters, adjusted to three frequency points. By 
applying the grey-body diagnostic to the ALMA band 7 
amplitude and slope, and including band 9, we find that 
the crescent is optically thick in band 7. Fig. [8] shows 
that grey-body line-of-sight models provide a fit to the 
observations, despite the use of a constant emissivity law 
fis(x) = 1.5. We used the CISscale flux-density scale; 
SEDscale produces similar maps but the low band 7 flux 
density results in lower optical depths. We also applied 
the grey-body diagnostic to the ATCA/ALMA multifre¬ 
quency data, as summarised in Fig. [9j The temperature 
map is essentially fixed by the band 9 spectral point. We 
note a minimum in temperature towards 0.5 h, which 
is also seen when using the ALMA data alone. The 
similarities between both temperature maps in CISscale 
is surprisingly good, considering the differences in uv- 
coverage, and the different treatments for f3. 

The uncertainties on the opacity index Ps(x) were es¬ 
timated by varying the ATCA intensities by Tier, where 
a is the scatter in the dirty map of image synthesis resid¬ 
uals. The resulting uncertainty map is shown in Fig. [9| 
The peak at 34 GHz, under the 11.5h clump, corresponds 
to /3s = 1.18 ±0.04 (in the SEDscale , similar significance 
values are obtained in the CISscale ), whereas the aver¬ 
age value outside the regions of low fts is 1.50, with a 
rms dispersion of 0.07 and typical values varying from 

I. 41 ±0.05 at lh, up to 1.63 ±0.08 at 9.5h. Compared to 
the average value of 1.50, the minimum in p s under the 

II. 5h clump is significant at ^7cr. 


F. SYNTHETIC DUST TRAP MODEL 


We constructed a synthetic model of a protoplanetary 
disk with a dust trap to make predictions comparable 
to the ALMA and ATCA observations of HD 142527. 
This synthetic model matches the observed SED (see be¬ 
low) , and is a perfected version of that whic h acc ounts for 
the scattered-light shadows (M arino et al.|2015] ), with a 
tilted inner disk. This model allows us to 1- assess if the 
observations can be interpreted with steady-state dust 
trapping formulae, given the tilted-inner disk, and 2- es¬ 
timate the biases involved in the uniform-slab grey-body 
model. 

The outer disk is where the dust trap is observed, so 
we describe it in more detail. It goes from 115 au to 
300 au with a rounded disk wall between 115 to 140 au, 
and is composed of 5.0xlO -5 M 0 in amorphous carbon 
grains and of 2.6 xlO -3 M 0 in silicate grains mixed with 
ices (using the Maxwell-Garnett rule). To simulate the 
dust trap effect on differen t grain sizes we used the ra¬ 
dial trap model from Pinilla et al. (|2012) cou pled with 
the p r escrip tions for azimuthal trapping from IBirnstiel 
et al. (2013). We first define an axisimmetric gas distri¬ 
bution with abounded disk wall, as previously described 


Marino et al. (2015), but with adjusted parameters 
7 = 7. w = 0. 13 and r r = 155.0 a u. We compute a criti¬ 
cal grain size (Pi nilla et al.||2012 ) a c at 133 au for which 


bigger particles would be radially trapped in the pressure 
maximum, while smaller particles would follow the gas 
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Figure 8. Spectral properties of the sub-mm continuum, obtained by comparing the intra-band ALMA band 7 spectral 
information with ALMA band 9, and with flux-density scale labelled CISscale. Upper panel: a optical depth map at 345 GHz, 
T 0 (x) b: flat emissivity law, taken here as f3 = 1.5 - this plot allows to better locate the asterisks that indicate the locations 
of the example spectra c: temperature map T s {x). Lower panel: example SEDs extracted at selected positions, indicated by 
asterisks in the upper panel. The y -axis shows specific intensity normalized by a v 1 spectrum, in arbitrary units. 


density background. Using a total gas mass of 0.25 M 0 
and an a t viscous parameter of 10 -2 we obtain a critical 
grain size of 0.7 cm at 133 AU. The Stokes number for 
0.7 cm grains, at peak gas densities of 2.6 10 10 cm -3 and 
20 K, is ^0.02 and close to a t . 

The above radial distribution is modulated in azimuth, 
using previ ously described formulae and symbols (Birn- 


stiel et ah 


20l3| (here we use R s = 148 AU, and" an 
azimuthal contrast in surface density of c = 20.0). The 
scale-height H is given by 


= 2a0 (iSTau) 


1.17 


(Fl) 


For the gas volume density we used the standard vertical 
Gaussian distribution. The small grains, i.e. amorphous 
carbon and silicate grains smaller than a c , follow the 
same density field as the gas. The azimuthal distribu¬ 
tion of big grains is calculated according to their Stokes 
number, a parameter to accoun t for the gas t urbulence a t 
and the gas background (Birnstiel et al.j 2013]). For their 
radial distribution we use a Gaussian parameterisation 


to account for the spread caused by turbulence: 

(r - Rs ) 2 


Pd,a{r > <l),z = 0) = Pd,a{4>) exp 


2<t 2 


> (F2) 


with a (a) = H r ^ ■ In order to compute the density 

as a function of height z, we use the same scale height 
as the gas. We investigated the implementation of grain 
settling, fol lowin g a parameterisation for the dust scale 
height ( Armitage|2010 ). However, we found that in order 
to reproduce the projection effect at ^lh, settling had to 
be very inefficient. Finally the distribution of big grains 
is normalized to the global grain size distribution and 

total dust mass in big g rains. _ 

We used RADMC3D ( Dullemond et al.|2015 ) to trans¬ 
fer the stellar radiation through the synthetic disk and 
predict emergent intensity maps, under the assumption 
of passive heating. We assumed a standard power-law 
distribution in grain sizes, with exponent —3.5. The 
outer disk is composed of three different dust species rep¬ 
resenting a range in sizes described by a m i n and a max : 










































Figure 9. Optical depth maps derived from the ATCA/ALMA multi-frequency data. The emissivity minimum at 11. 5h indicates 
the presence of larger grains. The upper panel corresponds to the absolute flux density scale labelled CISscale in Sec.[D] with 
a high flux density in band 7 (see Sec.|D|, while the lower panel corresponds to the flux scale in our preferred RADMC3D 
model, labelled SEDscale (see Sec.tDl). a: optical depth map at the reference frequency of 345 GHz. b: line of sight emissivity 
index map /3 s (x), with ATCA specific intensity contours in red. c: root-mean-square uncertainties on the emissivity index map. 
d: line of sight temperature, T s (x). 


• 5.0xlO -5 Mq of amorphous carbon grains with 
radii 1 < a < 100 /im with o ptical constants taken 
from Li & Greenberg (1997) 


• 2.5x10 3 M 0 in icy silicate grains, with radii of 
lOOjimi < a < a c and o ptic al constants from 


ning & Mutschke (1997) and Li & Greenberg 


Hen- 

998] ) 


• 1.0x10 4 M 0 in icy silicate grains, with radii o f 
a r < a < 1cm and optic al constants from Henning 


& Mutschke (1997) and Li & Greenberg ( 1998| ) 


As summarised in Fig.[l0j the predicted integrated flux 
densities fall close to thedoserved dereddened SED, ex- 


et al. 2011* 

Perez et al. 

2015b, 

archive. Fig. [1]J shows t 

hat th< 


1 


and from the herscke . 
re model also reproduces 


version of Fig. 2 in Marino et al. 2015). We switched-off 


scattering for the purpose 5T calculating images at sub- 
mm frequencies and lower. The reason for this is that 
the scattering cross sections calculated using Mie theory 
are probably overestimated, given that the actual grains 
are probably irregular in shape, or even porous, and so 
far from the Mie spheres. The emergent intensity maps 
in Fig. [T2| show that the 34 GHz emission exhibits a com¬ 
pact clump at the location of the dust trap, and a more 


extended crescent that follows the sub-mm morphology. 
In order to simulate the ATCA response, we sampled 
the i^-coverage of the actual ATCA observations, but 
on the radiative transfer predictions, and then averaged 
100 uvmem reconstructions with different realizations of 
noise. This exercise produced sky maps that are remark¬ 
ably similar to the ATCA-filtered observations. We also 
note that any flux loss is negligible. 

G. UNIFORM-SLAB DIAGNOSTICS ON THE 
SYNTHETIC DUST TRAP. 

We also applied the line-of-sight diagnostic based on 
grey-body models to the synthetic emergent intensities. 
As summarised in Fig. [l3j the basic features of the grey- 
body diagnostics approximate those seen in the data. 
There is a mild extension of the optical depth map 
into the northern ansa, which could be emphasized with 
higher optical depths. The center of the dust trap coin¬ 
cides with a minimum in /3 s- With this synthetic dust 
trap we can also test for biases in the grey-body diagnos¬ 
tics. The input optical depth map, 

/>00 

r(345 GHz)= / dspg(s) 

Jo 

/ da nj(a , g)^(a, 345GHz), (Gl) 
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Figure 10. Observed and model SED for HD 142527. 



Figure 11. Consistency of the synthetic dust trap model with 
the o bserved scatterd l i ght sh adows, a: NACO-PDI H-band image 
from Avenhaus et al. (2014) b: Radiative transfer prediction in 
H-band from tiie dust trap model. 


is remarkably consistent with the grey-body values - at 
native resolutions. We can also compare with the grey- 
body temperature diagnostic, T s (x), with the line-of- 
sight average temperature, as reported by RADMC3D 
after the thermal solution has converged, 


pOO 

(T)= dsp g {s) 

Jo 

J da^^T(a, s)rii(a, s) 


^(a, 345GHz) 
r(345GHz 


. (G2) 


Fig. [T3| shows that the two temperature diagnostics are 
very nearly equal (to within 10%). Limiting the integral 
to the depth s corresponding to r(345 GHz) = 1 resulted 
in essentially identical maps. 

Facilities: ATCA, ALMA. 
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Figure 12. Predicted emergent continua at 34 GHz, 345 GHz and 700 GHz. Upper panel: raw RADMC3D prediction. Bottom 
panel: predictions filtered by ATCA+uvmem response, with restoration in natural weights. Note the 11.5h clump of cm-sized grains, which 
stands out at ATCA frequencies, but is absent at ALMA frequencies. 
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Figure 13. Grey-body ’line-of-sight’ diagnostics calculated on the radiative transfer predictions, and comparison with input 
optical depth and temperatures. A mask has been applied to compute line-of-sight diagnostics at a fraction of 1/5 times peak 
intensity, a), b), c) : grey-body diagnostics calculated after filtering by the ATCA+uvmem response - a) shows the optical depth 
map at the reference frequency of 345 GHz, b) shows the line of sight emissivity index map /3 s (x), and c) shows line of sight 
temperature, T s (x). d), e), f): same quantities, as calculated on the radiative transfer predictions, g), h): optical depth and 
average temperature maps used as input for the ray-tracing. See text for the definition of average temperature. 






















































